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Обсуждаются вопросы создания математических моделей биологических объектов, которые необходимы при проведении 
электроимпедансных томографических исследований. Предложен критерий оптимизации размеров конечных элементов моде- 
ли и проведены оценки максимального количества информации, получаемого при различных параметрах аппаратуры для 
электроимпедансных измерений. На примере построения модели головы человека рассмотрены приемы моделирования био- 
логических объектов. 


Электроимпедансная томография (ЭИТ ) это 
метод реконструкции распределения проводимос- 
ти внутри объекта на основе результатов электри- 
ческих измерений на поверхности. Если каждый 
электрод имеет площадь 5, м 2 , контактное сопро- 
тивление Д,, Ом и через него проходит ток I, А, то 
полная модель измерения в электроимпедансной 
томографии выглядит следующим образом: 

Ѵ(сгѴф) = 0; 

Ѵ А =(р + К сХ +а^-\ 
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где ст - удельная проводимость, См/м, Ѵ е - напря- 
жение на электроде относительно нуля, вольт, по- 
тенциал электрического поля. Нижние индексы от- 
носятся к номеру электрода в паре. Для измеритель- 
ных электродов модель будет такой же, но ток I бу- 
дет равен нулю, соответственно, подинтегральные 
выражения в (1) тоже будут равны нулю. Выраже- 
ние (1) - это постановка прямой задачи электроим- 
педансной томографии - найти распределение по- 
тенциала электрического поля при заданном расп- 
ределении проводимости и расположении электро- 
дов, подводящих ток. Несмотря на то, что фунда- 
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ментальные соотношения, описывающие поведе- 
ние электромагнитного поля в биологических объ- 
ектах, в некоторых случаях имеют точное, аналити- 
ческое решение, для практической реализации ме- 
тода ЭИТ необходимы численные методы нахожде- 
ния решения прямой и обратной задач. В основе 
всех численных методов лежит приближение, ап- 
проксимация, пространственно непрерывного 
распределения проводимости, его дискретным 
представлением в пределах границ изучаемого объ- 
екта. Удельное сопротивление биологического объ- 
екта на макроскопическом уровне распределено 
непрерывно в зависимости от координат, но мы не 
можем точно описать его каким-либо математичес- 
ким выражением и потому вынуждены прибегнуть 
к определенному упрощению. Для этого необходи- 
мо создать модель исследуемой области. Основная 
цель создания математической модели состоит в 
том, чтобы максимально упростить реконструкцию 
изображения при минимальной потере точности. 

Необходимость и возможность создания приб- 
лиженной модели следует из того, что реально про- 
водимые в электроимпедансной томографии изме- 
рения дискретны как во времени, так и в простран- 
стве. В самом деле, для наблюдения любой величи- 
ны с помощью аналого-цифровых преобразовате- 
лей (АЦП) всегда требуется определенный, пусть да- 
же очень малый, интервал времени. Кроме того, для 
получения точечных статистических оценок изучае- 
мого процесса чаще всего требуется, чтобы за время 
измерения не произошло сколько-нибудь сущест- 
венного изменения электрических параметров, поэ- 
тому интервал времени регистрации импеданса би- 
ологического объекта обычно составляет от 10 5 до 
10 1 с. При уменьшении времени измерения появля- 
ются достаточно серьезные технические проблемы. 
Современные АЦП позволяют достичь времени од- 
ного измерения не более долей мкс при отношении 
сигнал/шум до 80 дБ, поэтому они все чаще приме- 
няются для так называемой "цифровой квадратур- 
ной демодуляции" с последующей цифровой фильт- 
рацией [1]. Однако уменьшение времени измерения 
уменьшает отношение сигнал/шум и увеличивает 
погрешности измерения. Конечно, существует воз- 
можность использования наносекундных АЦП, но 
применять их в электроимпедансной томографии 
нецелесообразно по двум причинам. Во-первых, 
резко возрастает стоимость аппаратуры. Кроме того, 
скорость протекания наиболее важных физиологи- 
ческих процессов, влияющих на электрические 
свойства биологических тканей, относительно неве- 
лика верхняя граница спектра реографического сиг- 
нала не превосходит 10.. .20 Гц. Поэтому даже в 
электроимпедансных системах реального времени с 
частотой регистрации 20-30 кадров в секунду при 16 
электродах минимальное время измерения состав- 
ляет десятки и сотни мкс и уменьшать его нет ника- 
кой практической необходимости. 

Увеличение времени измерения сверх долей се- 
кунды приводит к необходимости обеспечивать 
постоянные условия регистрации в течение про- 


должительного времени (несколько десятков с), 
что существенно затрудняет проведение исследо- 
ваний іп ѵіѵо и дает возможность изучать только 
статические или очень медленные процессы. 

Кроме дискретности во времени, существует и 
дискретность в пространстве, поскольку для прове- 
дения неинвазивных измерений доступна только 
поверхность биологического объекта и только в ог- 
раниченном числе точек. Как правило, в существу- 
ющих системах применяется от 16 до 256 электро- 
дов. Но, сколько бы электродов не использовалось, 
все равно сохраняется пространственная дискрети- 
зация исследуемой области с большим или мень- 
шим межэлектродным расстоянием. Таким обра- 
зом, в реальных условиях нам доступны дискретные 
измерения импеданса через определенные проме- 
жутки времени в некоторых точках поверхности би- 
ологического объекта. Интуитивно понятно, что 
погрешность дискретизации (как временной, так, и, 
в особенности, пространственной) имеет достаточ- 
но большую величину и нет смысла строить матема- 
тическую модель с более низкой погрешностью, чем 
погрешность дискретизации. В данной работе мы 
предлагаем новый подход, на основании которого 
можно строить математические модели для электро- 
импедансной томографии и оценивать их качество. 

В основу нашей идеи положен информацион- 
ный подход. Пусть имеется система с М электродов 
и отношением сигнал/шум при измерении импе- 
данса А, что соответствует Г~А/6-разрядному 
АЦП. Количество пространственных линейноне- 
зависимых измерений для получения одного кадра 
изображения составляет К=М(М-Н)/ 2. Постоянная 
И может принимать счетные значения от 1 до 3 в за- 
висимости от выбранного способа коммутации 
электродов. Оценим количество информации, ко- 
торое в пределе можно получить с помощью дан- 
ной системы. Одно линейно-независимое измере- 
ние создает I бит информации. Количество таких 
измерений в кадре - К. Используя комбинаторное 
определение количества информации по Хартли 
[2], получим очевидное соотношение: 

У =ЬК1о% 2 2=ЬК. (2) 

Если используется двоичное представление из- 
меренной величины, то количество возможных зна- 
чений одного измерения составляет 2 1 . Эти комби- 
нации образуют алфавит, состоящий из 7 1 символов. 
Соответственно, ^измерений приведут к появлению 
2“ комбинаций. После логарифмирования получаем 
(2). Рассмотрим конкретный пример. Пусть имеется 
система с 16 электродами и отношение сигнал/шум 
составляет 70 дБ. Измерение осуществляется при по- 
даче тока между двумя соседними электродами па- 
рой двух других электродов, на токовых электродах 
измерения не проводятся. В этом случае Н= 3, М= 16, 
Х=104, 1=12. Эти параметры характерны для извест- 
ных томографических систем [3]. Количество ин- 
формации после проведения всех измерений для од- 
ного кадра будет равно 7=1248 бит. Предельная раз- 
решающая способность подобной измерительной 
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системы в точности определяется количеством ли- 
нейно-независимых измерений в кадре, но, в отли- 
чие от рентгеновской или магнито-резонансной то- 
мографии, размер элемента изображения нелинейно 
зависит как от проводя щ их свойств участка исследу- 
емого объекта, так и от измерительной схемы и рас- 
положения электродов. Уверенно можно говорить 
лишь об одном количество информации в рекон- 
струированном изображении не может быть больше 
количества информации, полученной в результате 
измерения. Например, если в итоговом изображении 
системы из нашего примера предполагается отобра- 
жать значения проводимости 256 градациями серого 
(8 бит на элемент изображения), то модель исследуе- 
мой области может содержать не более 156 элемен- 
тов. Еще одна особенность электроимпедансных из- 
мерений состоит в том, что напряженность электри- 
ческого поля максимальна на периферии и резко 
(пропорционально кубу расстояния от поверхности) 
убывает к центру исследуемой области. 



Рис. 1. Поверхность сферической модели головы 



Рис. 2. Зависимость напряжения в узлах модели от расстоя- 
ния до токового электрода 


Рассмотрим это на конкретном примере. Пусть 
имеется четырехслойная сфера (рис. 1) с удельной 
проводимостью слоев от 0,042 (кость) до 4,6 (ликвор) 
См/м и диаметром 0,15 м. На границе области выде- 
лено 17 подобластей, в которых размещаются "элект- 
роды”. Каждый электрод моделируется совокуп- 
ностью поверхностных треугольных элементов. Все- 
го модель содержит 5187 конечных элементов в фор- 
ме тетраэдра. На один из "электродов" подается ну- 
левой потенциал для обеспечения единственности 
решения, на диаметрально противоположный пода- 
ется ток силой 1 мА. В результате решения прямой 
задачи вычисляются значения потенциала электри- 
ческого поля в узлах модели. Для придания физичес- 
кого смысла решению вместо потенциала нужно ис- 
пользовать разность потенциалов в узлах по отноше- 
нию к референтному уровню на нулевом электроде. 
На (рис. 2) приведена зависимость разности потен- 
циалов от расстояния до токового электрода. 

Вклад центральных элементов в результат изме- 
рения существенно ниже по сравнению с перифери- 
ческими, соответственно, они дают меньший вклад в 
общее количество информации, создаваемое при из- 
мерениях. Пространственное разрешение электро- 
импедансного томографа зависит от его способности 
различать малые изменения разности потенциалов 
за счет различий в проводимости отдельных участков 
исследуемой области. Из (рис. 2) интуитивно ясно, 
что наибольшие колебания разности потенциалов 
наблюдаются возле поверхности и убывают к центру. 
Для количественной оценки этого явления можно 
использовать так называемую функцию чувстви- 
тельности [4], которая сопоставляет каждой точке 
исследуемой области наблюдаемую на периферичес- 
ких электродах разность потенциалов при заданной 
конфигурации зондирующего тока. Эта функция в 
точности равна найденному в результате решения 
прямой задачи потенциалу электрического поля. Но 
для электроимпедансной томографии необходимо 
знать, как неоднородность удельной проводимости 
повлияет на результат измерений: 

8.1 = У^басрсрсІО.-, ( 3 ) 

О. 

где / - вариация результата измерений, ср* и ср - 
функции чувствительности соответственно токо- 
вого и потенциального электродов, да - вариация 
проводимости, сІО. - минимальный размер элемен- 
та исследуемой области (площадь или объем). 

Функция чувствительности каждого электрода 
определяется путем решения прямой задачи, при 
условии что на данный электрод подается ток. Оче- 
видно, что для получения электроимпедансного 
изображения нужно, как минимум, регистрировать 
<5/. Это дает нижнюю оценку диапазона измеряе- 
мых величин разности потенциалов между поверх- 
ностными электродами. Погрешность измерения е 
не должна быть больше $=да(р* срсІ ( Л. При задан- 
ной точности измерительного тракта и геометрии 
исследуемой области, величину 5 можно изменять 
двумя способами - изменением размеров конеч- 
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ных элементов и функции чувствительности. Уве- 
личения чувствительности можно добиться путем 
выбора оптимальной конфигурации зондирующе- 
го тока [5]. Однако, этот способ имеет существен- 
ные ограничения, связанные как с усложнением 
аппаратуры, так и с геометрией исследуемой облас- 
ти. Вообще говоря, оптимальная конфигурация 
зондирующего тока единственна для заданного 
распределения источников тока и удельной прово- 
димости, но для ее нахождения необходимо знать 
распределение удельной проводимости. Таким об- 
разом, для полноценной реализации этого способа 
повышения чувствительности необходимо увязы- 
вать между собой алгоритм реконструкции изобра- 
жения и алгоритм регистрации, но это весьма 
сложная задача и на сегодняшний день оптималь- 
ная конфигурация зондирующего тока рассчитана 
только для простых областей (окружности) при ус- 
ловии, что вариации удельной проводимости неве- 
лики. Кроме того, при такой конфигурации зонди- 
рующего тока возникают чисто технические проб- 
лемы уменьшения контактного сопротивления 
электродов и обеспечения стабильности условий в 
месте контакта электрод-кожа. 


Из (3) следует, что 5” можно увеличить и прос- 
тым изменением размера элемента исследуемой 
области сЮ.. Так как всегда можно соблюсти усло- 
вие 0<аЪ$(ср* ср(10.)<\ , то можно найти величину: 


р = ^-г1°в(г); где 

П 


г = аЪ$((р срсЮ). 


(4) 


По аналогии с информационной энтропией, ве- 
личину р тоже будем называть энтропией. Наиболь- 
шее значение р будет соответствовать максимально 
однородному пространственному распределению 
чувствительности внутри исследуемой области. Вы- 
ражение (4) дает возможность найти наилучшее 
распределение размеров элементов путем нахожде- 
ния экстремума. В том случае, если погрешности из- 
мерительного тракта обусловлены, в основном, вли- 
янием независимых источников шума, такой выбор 
размеров элементов модели области интереса поз- 
воляет сузить динамический диапазон измеряемых 
величин (либо увеличить точность измерения при 
заданном диапазоне измерительного тракта) и, тем 
самым, наиболее полно использовать возможности 
аппаратной части. Кроме того, ограничивается 
сложность математической модели и существенно 
уменьшаются вычислительные затраты на рекон- 
струкцию изображения. В целом предлагаемый на- 
ми подход позволяет оптимальным образом согла- 
совать свойства измерительного тракта и разреше- 
ние реконструированного изображения. Единствен- 
ный очевидный недостаток связан с необходи- 
мостью решения нелинейной задачи оптимизации 
при нахождении размеров элементов, но для совре- 
менных вычислительных средств это вполне по си- 
лам. Существует большое количество программ для 
реализации метода "максимальной энтропии" в дру- 
гих областях, в частности, свободно распространяе- 


мый пакет Кещііашаііоп Тоо1$ (Нашей Р.С., 
ййр://\ѵ\ѵ\ѵ.іптт.йІи.йк/~рс1т). Таким образом, при 
реконструкции распределения проводимости мож- 
но существенно увеличить допустимую погреш- 
ность нахождения проводимости центральных эле- 
ментов. Это эквивалентно уменьшению количества 
информации, приходящегося на центральные об- 
ласти. В результате открывается возможность адап- 
тивного изменения размеров элементов изображе- 
ния - центральные элементы можно сделать больше 
(уменьшить пространственное разрешение), а пери- 
ферические -меньше (увеличить разрешение). Ин- 
терес к адаптивным методам построения конечно- 
элементной модели (КЭМ) проявился недавно и в 
настоящее время это весьма многообещающие нап- 
равление повышения разрешающей способности 
электроимпедансной томографии [6]. 

Основой для применения адаптивных методов 
модификации КЭМ служит исходная модель, пост- 
роенная, может быть, с избыточной точностью. Во- 
обще говоря, существует несколько способов приб- 
лиженного описания пространственного распреде- 
ления электрических свойств биообъектов - с по- 
мощью метода конечных элементов, метода конеч- 
ных разностей, метода граничных элементов, раз- 
личных полуаналитических способов и т.д. Не вда- 
ваясь в тонкости реализации каждого из известных 
подходов, отметим лишь, что наибольшее распро- 
странение в электроимпедансной томографии по- 
лучил метод конечных элементов, именно он обес- 
печивает достаточную для практических целей точ- 
ность моделирования и в то же время приемлемую 
вычислительную сложность. 

Последовательность создания конечно-элеме- 
нтной модели предметной области состоит из нес- 
кольких этапов. Вначале необходимо создать кон- 
турную модель границы моделируемого объекта. 
Это можно сделать различными способами. Наибо- 
лее простой состоит в создании виртуальной моде- 
ли в одной из систем трехмерного моделирования. 

Для проведения предварительных исследований 
мы использовали модель головы, входящую в состав 
примеров к системе автоматизированного проекти- 
рования СЮ (демонстрационная версия и описание 
доступны по адресу Шр://§іс1.сітпе.ирс.е$). Вначале 
модель представляет собой набор линий и поверх- 
ностей, описывающих внешний контур органа или 
участка тела. Затем осуществляется разбиение пове- 
рхности на треугольники (триангуляция) и получает- 
ся набор конечных элементов, с помощью которых 
аппроксимируется внешняя граница моделируемой 
области. После создания границы можно сформиро- 
вать разбиение внутренней области. Для этого опре- 
деляется объем, ограниченный поверхностью объек- 
та и генерируется разбиение на конечные элементы, 
имеющие, например, форму тетраэдра (рис. 3). За- 
тем необходимо присвоить каждому конечному эле- 
менту значение проводимости. Это наиболее трудо- 
емкая операция, особенно при большом количестве 
элементов. В большинстве доступных систем авто- 
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матизированного проектирования задание свойств 
конечных элементов можно производить в пакетном 
режиме, либо вручную, используя специальную мас- 
ку, накладываемую на модель. В любом случае, для 
создания адекватной модели может потребоваться 
очень много времени и ручной работы. В наших ис- 
следованиях нам потребовалось около 200 ч рабоче- 
го времени для создания близкой к реальности моде- 
ли головы. В качестве эталона были использованы 
сегментированные томографические изображения, 
на которых были выделены отдельные ткани. 



Рис. 3 . Разбиение трехмерной модели головы на конечные 
элементы 

Качество полученной модели определяется ее 
способностью моделировать исследуемой процесс 
с заданной степенью точности. По сравнению с 
традиционно используемыми моделями в виде со- 
вокупности вложенных сфер, точность моделиро- 
вания распределения электрического потенциала 
увеличивается на 15... 20 % при условии масштаби- 
рования модели в соответствии с реальными разме- 
рами головы. Но нам встречались пациенты, фор- 
ма головы которых значительно отличалась от мо- 
дельной и погрешность моделирования в этом слу- 
чае была сравнима со сферической моделью. Не- 
достаток данного подхода проявляется, прежде 
всего, в его трудоемкости и высокой стоимости 
систем автоматизированного проектирования, без 
которых обойтись весьма сложно. 

Другой подход был реализован [7] с помощью 
системы ШЕА8. За основу в этой работе были взя- 
ты томографические магнито-резонансные изоб- 
ражения (можно также использовать и рентгеновс- 
кие). Толщина среза должна составлять 3...5 мм, 
максимум 10 мм, в противном случае неизбежен 
эффект "ступенчатой" границы. Срезы должны ох- 
ватывать всю моделируемую область с небольшим 
(2-3 дополнительных среза) запасом. Каждый срез 
подвергается дополнительной обработке с целью 
точного выявления границы всего объекта и гра- 
ниц анатомических структур, отличающихся про- 
водимостью. В простейшем случае это выполняет- 
ся с помощью фильтра высоких частот, но для кор- 


ректной сегментации требуются значительно более 
сложные алгоритмы. Вопросы сегментации изоб- 
ражений выходят далеко за рамки нашей статьи и 
мы не будем их рассматривать, тем более, что суще- 
ствует огромное количество как печатных работ, 
так и доступных программ с открытыми исходны- 
ми текстами для выделения контуров и сегмента- 
ции изображений (например, комплекс №да§е8, 
созданный ВсгпИагсІ Оеі§ег и Іеап-ОапіеІ Воікзоппаі 
Соругцфй (С) 1993-2001 Ьу ПЫКІА Зоріііа Апііроііх). 
В состав №іа§е$ входит набор различных программ, 
в том числе и для сегментации, распространяется 
для некоммерческого использования в виде исход- 
ных текстов. В результате сегментации получается 
набор контуров в виде отрезков прямых линий или 
сплайнов более высоких порядков. Непосред- 
ственно использовать полученные контуры неце- 
лесообразно, поскольку они содержат слишком 
много мелких деталей, которые могут существенно 
ухудшить качество реконструкции поверхности би- 
ообъекта. Для упрощения контура также существу- 
ют специализированные программы (в Миа§е§ та- 
кой инструмент имеется). Затем из совокупности 
контуров реконструируется поверхность моделиру- 
емого объекта. Единственное требование к создан- 
ной поверхностиона должна быть замкнутой и не 
содержать разрывов. Поверхность создается в виде 
триангуляции (совокупности треугольников), на 
основе которых можно затем осуществить разбие- 
ние внутренней области на конечные элементы. 
Эта задача решается с помощью так называемого 
генератора сеток (те$Ь §епегаІог). Его основная за- 
дача - разделить внутренний объем, ограниченный 
поверхностью, на конечные элементы (например, 
тетраэдры) с заданными свойствами. В своей рабо- 
те мы успешно применяли программу ТеіСеп (ав- 
тор Нап§ 8і, Е-таі1: 8І@\ѵіа8-Ьег1іп.сіе), задавая ис- 
ходные данные в виде многоугольников (роіу файл 
в терминологии ТеіСеп). Сегментация позволяет 
найти отдельные области и для каждой из них 
построить КЭМ. В результате получается набор от- 
дельных моделей, которые нужно свести воедино. 
Это можно сделать либо с помощью специализиро- 
ванной программы (например, 8е§Ма1с1г, ОЫо 
8іаІе Шіѵегеііу), либо непосредственно в среде ви- 
зуального проектирования. В работе К.Н. ВаѵТогсІ 
еі.аі. использована универсальная среда ГОЕА8, но 
можно использовать и более доступные и распро- 
страненные средства, такие, как описанный выше 
пакет ОШ. После объединения всех частей в одну 
получается законченная конечно-элементная мо- 
дель исследуемой области. Значения проводимости 
каждому элементу нужно присвоить после сегмен- 
тации до сведения отдельных частей воедино, пос- 
кольку в пределах каждой части проводимость при- 
нимает одно и то же значение. 

Вне зависимости от способа построения конеч- 
но элементной модели на завершающем этапе нуж- 
но выделить точки фиксации электродов в соответ- 
ствии с определенной схемой. Мы использовали 
общепринятую схему "10/20", широко применяе- 
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мую в традиционной электроэнцефалографии. На- 
иболее сложный для моделирования вопрос - это 
учет погрешности размещения электродов во время 
проведения реальных исследований по сравнению 
с моделью. Ошибки в размещении электродов при- 
водят к существенным геометрическим искажени- 
ям изображения и оказывают не меньшее влияние 
на конечный результат, чем погрешность измере- 
ния напряжения на электродах. Разумный компро- 
мисс между сложностью проведения исследования 
и ошибками реконструкции можно достичь, изме- 
ряя расстояние от кончика носа до затылочного 
бугра обследуемого, а также расстояние между слу- 
ховыми проходами. КЭМ легко масштабируется в 
соответствии с заданными размерами и погреш- 
ность можно уменьшить на несколько процентов. В 
действительности все зависит от формы головы и 
погрешности установки электродов. Какие-либо 
теоретические оценки сделать практически невоз- 
можно. Во время создания КЭС мы пытались изме- 
рить погрешность аппроксимации головы путем 
сравнения магнито-резонансной томограммы и 
среза КЭМ. Для этого мы накладывали срезы друг 
на друга и измеряли площадь несовпадающих час- 
тей. Измерения проводились с помощью ориги- 
нальной программы для обработки изображения, 
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созданной на кафедре медицинской и биологичес- 
кой кибернетики Сибирского государственного ме- 
дицинского университета. Полученная величина 
составляла от 0,1 до 10 % от площади среза. 

Таким образом, в результате создания конечно 
элементной модели получается приближенное (в 
пределах определенной погрешности) математичес- 
кое описание исследуемой области. Модель состоит 
из точек (узлов), в которых соединяются линии, об- 
разующие конечные элементы. Описание модели 
содержит, как правило, три основных компонента - 
перечень всех узлов с указанием их пространствен- 
ных координат, перечень конечных элементов с ука- 
занием входящих в каждый элемент узлов и значе- 
ния проводимости для каждого элемента. 

Создание модели - это один из самых трудоем- 
ких этапов в разработке системы для электроимпеда- 
нсной томографии. Для предварительных исследо- 
ваний возможно использование простых описаний, 
которые моделируют исследуемый объект одной или 
несколькими геометрическими фигурами (сфера, 
эллипсоид, куб, цилиндр и Т.Д.). Такие модели с по- 
мощью специализированных программ можно соз- 
давать очень быстро и они могут вполне адекватно 
описывать поверхность отдельных участков тела. 
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